Method for measuring the entanglement of formation for arbitrary-dimensional pure 

states 
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Entanglement of formation is an important measure of quantum entanglement. We present an 
experimental way to measure the entanglement of formation for arbitrary dimensional pure states. 
The measurement only evolves local quantum mechanical observables. 
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Quantum entangled states have become the most im- 
portant physical resources in quantum communication, 
information processing and quantum computing. One 
of the most difficult and fundamental problems in en- 
tanglement theory is to quantify the quantum entangle- 
ment. A number of entanglement measures such as the 
entanglement of formation and distillation [JJ, Q , negativ- 
ity [|[ and concurrence [H have been proposed. Among 
these entanglement measures, the entanglement of forma- 
tion, which quantifies the required minimally physical re- 
sources to prepare a quantum state, plays important roles 
in quantum phase transition for various interacting quan- 
tum many-body systems [f| and may significantly affect 
macroscopic properties of solids @. Thus the quantita- 
tive evaluation of entanglement of formation is of great 
significance. 

Comparing with the concurrence, entanglement of for- 
mation is more difficult to deal with, and less results have 
been derived. The experimental measurement of concur- 
rence has been proposed in 0] for pure two-qubit sys- 
tems by using two copies of the unknown quantum state. 
In [j| the authors presented an approach of measuring 
concurrence for arbitrary dimensional pure multipartite 
systems in terms of only one copy of the unknown quan- 
tum state. However, due to the complicated expression 
of entanglement of formation, there is no experimental 
way yet to determine the entanglement of formation for 
an unknown quantum pure state, except for the case of 
two-qubit systems for which the concurrence and entan- 
glement of formation have a simple monotonic relations 



In this brief report we present an experimental de- 
termination of the entanglement of formation for arbi- 
trary dimensional pure quantum states. The measure- 
ment only evolves local quantum mechanical observables 
and the entanglement of formation can be obtained ac- 
cording to the mean values of these observables. 

The entanglement of formation is defined for bipartite 
systems. Let Ha and Hb be m and n (m < n) dimen- 
sional complex Hilbert spaces with orthonormal basis \i), 
i = l,...,m, and i = l,...,n respectively. A pure 



quantum state on Ha <8> Hb is generally of the form, 

m n 

H) = Y Y ai J ^ ' aij ec (1) 

i=l 3 = 1 



with normalization 



Y Y °« a « = 1 

i=l j=l 



(2) 



The entanglement of formation of is defined as the 
partial entropy with respect to the subsystems 

E{m = -Tr (p A log 2 p A ) = -Tr (p B log 2 p B ) , (3) 

where p A (resp. p B ) is the reduced density matrix ob- 
tained by tracing \ip)(ip\ over the space Hb (resp. Ha)- 
This definition can be extended to mixed states p by the 
convex roof, 



%)e min Vp,£(|i)), 
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(4) 



where the minimization goes over all possible ensemble 
realizations of p, 



P : 
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(5) 



A bipartite quantum state \ip) can be written in the 

m 

Schmidt form \ip) = J2 |m)|«b)> K > 0, = h 

under suitable basis \%a) € Ha and \%b) € Hb- K, i — 
I,...,m, are also the eigenvalues of p A . E(\tfj)) can be 
further expressed as 



E(m = S(p A ) = -J2x i log Aj. 



(6) 



For two-qubit case, m = n = 2, \ip) = an|00) + 
ai2|01)+a 2 i|10)+a 22 |ll)Jaii| 2 + |ai 2 | 2 + |a 2 i| 2 -|-|a 2 2 12 
I. © can be written as Q, 



(7) 
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where h(x) 



-x log 



(1 



ilog 2 (l 
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2|ana22 — ai2Q2i| is the concurrence. In this special case 
E{ | ip)) is just a monotonically increasing function of the 
concurrence C. However for m > 3, there is no such rela- 
tions like (O between the entanglement of formation and 
concurrence in general. Since for the case to = 2, due to 
the normalization condition, Ai + A2 = 1, only one free 
parameter is left in the formula ([5]). For general high 
dimensional case, E(\i/})) depends on more free param- 
eters. Nevertheless, if p A has only two non-zero eigen- 
values (each of which may be degenerate) , the maximal 
non-zero diagonal determinant D of p A is a generalized 
concurrence, namely, the corresponding entanglement of 
formation is again a monotonically increasing function of 
D M. The construction of such kind of states is presented 
in [lOj ]. In [ll|, the results are generalized to more gen- 
eral case: relations like ([7]) holds for states with p A having 
more non-zero eigenvalues such that all these eigenvalues 
are functions of two independent parameters. 

To measure the quantity ^ experimentally, we first 
rewrite the expression ^ according to the entanglement 
of formation of some "two-qubit" states. Let L A and 
Lp be the generators of special unitary groups SO{m) 
and SO(n), with the m(m — l)/2 generators L A given 
by {|«)0'l ~ |j><*|} 3 I < i < j < m, and the n(n - l)/2 
generators Lp given by — 10(^1}, t < k < I < n, 

respectively. The matrix operators L A (resp. Lf ) have 
rn— 2 (resp. n~ 2) rows and m — 2 (resp. n — 2) columns 
that are identically zero. 

Let p — \il)){ip\ be the density matrix with respect to 
the pure state We define 
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1,2, 



-i(n-l) 



(8) 



and 



where a 

||X|| = y/Tr(XX^) is the trace norm of matrix X. As 
the matrix L A ® La has mn — 4 rows and ran — 4 columns 
that are identically zero, p a p has at most 4x4= 16 
nonzero elements and is called "two-qubit" state. p a p is 
still a normalized pure state. 

Theorem 1 For any m®n (to < n) pure quantum state 



E(p a p) + \og{C a p) 

C a /3 



, (9) 



where C afj = l/Tr{L A <g> Lj\^) {i>\{L A y ® (Lf )*}. 

Proof. To calculate E{p a p) we denote = |a)(6| — 
1 6) (a. I and = |c)(d| — \d)(c\ for convenience, where 
1 < a < b < m and 1 < c < d < to. Set 

^ = ia®^flV')^l(^) t ®(£f) t . (10) 

It is direct to verify that 



Pa/3 



)af3{ 



(11) 



where | ip 
X a Sac\bd). 



Ct0 



XbSbd\ac) - X b 5 bc \ad) 



We now compute the eigenvalues of pA 



Xa$ad\bc) + 
= Tr B{p' a0 ) 



according to the values of a, b, c and d: 

i) . a ^ b ^ c ^ d. We have \ip) a p — 0. 

ii) . b > a = c < d and b 7^ d. We get \ip) a p = y/\^\bd). 
The eigenvalue of p A p corresponding to this case is A a - 
Asa = c can be chosen to be 1, 2, • • • , m — 2, b and d have 
only to — fc and to — k — 1 (corresponding to a = c = k, 
fc = l,2,---,m — 2) kinds of choices. Altogether we have 
(to — k) (to — fc — 1) eigenvalues of to be Afc in this 
case, with k = 1,2, • • • ,m — 2. 

iii) . a < 6 = rf > c and a ^ c. We have |'0)q / 9 = 
y/\~b\ac). The eigenvalue of is A/,. In this case b = d 
can be 3, 4, • ■ • , to. Then a and c have only k — 1 and 
fc — 2 (corresponding to 6 = <i = A;, k = 3, 4, • • • , to) kinds 
of choices. Hence we have (k — l)(k — 2) eigenvalues of 
pA to be Afc in this case, k = 3, 4, • • • , to. 

iv) . 6 > a — c < d and b = d. We obtain \ip) a p = 
-\/Afc|ac) + \fX^\bd). The eigenvalues of are A Q and 
Af,, and a = c can be 1, 2, • • • , to — 1. Then b = d can 
be k + 1, fc + 2, • • ■ , to (corresponding to a = c = fc, fc = 



1) eigenvalues of to 



1, 2, • • • , to — 1). We have (to 
be Afc, k — 1, 2, • • • , to. 

v). a < b = c < d. We have IV'W = —\f\~b\ad). The 
eigenvalue of p A p is Ab. b = c can be 2,3, •• ■ , m — 1. 
Then a and <i have only k — 1 and m — k (corresponding 
to & = c = fc, fc = 2, 3, • • • , to — 1) kinds of choices. 

k) eigenvalues of pA to be Afc, 



We have (k — 1)(ti 
k = 2,3,--- ,m- 1. 

vi). c < d = a < 6. We have \^) a p — — \f\~2\bc). 
The eigenvalue of pA is A a - In this case a = d can 
be 2, 3, • • • , to — 1. c and b have only fe — 1 and m — fc 
(corresponding to 6 = c = fc, fc = 2, 3, • • • , to — 1) kinds 
of choices. Therefore we have (fc — l)(m — fc) eigenvalues 
of p A p to be Afc, with fc = 2, 3, • • • , m — 1. 

Let \ l a p stand for the eigenvalues of p A p. From the 
analysis of cases i)-vi) and formula ([6]) , we get 



(m - 1) 



^^A^log(A^) 

a/3 i=l 



(12) 



Since 



Pa/3 n 



(13) 



Tr{p' a[j } 

we have ^ ^L^Cq^ = 1 f° r an y a an d P- Therefore 

i=l 

i=l i=l 

= - log(C Q/3 ) - A^ log(A^). 
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That is 



^A^log(A^) 

i=l 



Ejpap) + \og{C a p) 



Substituting (iPtl) into (fT2")l . we obtain that 



(m — 1) 



l) 2 ^ 



E{p a p) + log(C Q/3 ) 



a/3 



a 



a/3 



(14) 



(15) 



which proves the theorem. □ 

The theorem shows that one can derive the entangle- 
ment of formation of a pure quantum state by measur- 
ing the values of the entanglement of formation of all 
the states p a p and values of C a p. Here if \ip) a p = 0, 
then C a p goes to infinity and this term does not con- 
tribute to the summation in (jH]). Hence the summation 
J2a/3 m © simply goes over all the terms such that 
Tr{Li® Lf|^|(^)t®(Lf)t}^0. 

With formula ©, we now show how to get the value 
of E{\tl))) experimentally by measuring the quantities on 
the right hand side of ©. 

The quantity C a p — l/Tr{p a/3 } can be determined by 

Tr{p' ap }. Since Tr{p a0 } = (ip\(L^L^ ® {Lf^Lf\f), 
one can obtain C a p by measuring the local Hermitian 
operator (L^)^L^ (g> (La )^£f associated with the state 

To measure E(p a p), we first note that although p a p are 
m®n bipartite quantum states, they are basically "two- 
qubit" ones. For given L a — — and Lp — 

\k)(l\ — \l)(k\, i ^ j, k 7^ I, the non-zero elements of p a p 
are located at the i * (m — 1) + fcth, i * (m — 1) + ith, 
j * (m — 1) + fcth, and j * (m — 1) + Zth rows and the 
i * (m — 1) + fcth, i * (m — 1) + ith, j * (m — 1) + A:th, 
and j * (rn — 1) + /th columns. They constitute a 4 x 4 
matrix, 



'a/3 



( Pik,ik Pik,il Pik,jk Pik,jl 

Pil,ik Pil,il Piljk Piljl 

Pjk,ik Pjk,il Pjk.jk Pjk.jl 

V Pjl.ik Pjl.il Pjl.jk Pjl.jl 



Set a a p = v' a p/Tr{<j' al3 }. Obviously E(p a p) = E(a a p). 
But a a p are actually two-qubit pure states. According 
to the formula ([7]), E(a a p) is determined by the concur- 
rence C(a a p) = C(p a p). Therefore if we can measure the 
quantity C(p a p), we can obtain E(p a p). 



The quantity C(p a p) can be measured experimentally 
in terms of the method introduced in Q , with a few mod- 
ifications of the measurement operators. Corresponding 
to the case of L a = — and Lp = \k)(l\ — \l){k\, 
we define m x m matrix operators E s , s = 0,l,2,3, such 

that (Eo)pg = SpiSqi + Spjfiqj, {^l)pq — dpifiqj + SpjSqi, 
fi- J 2)pq — I&pi^qj I^pj^qi, \^s)pq — ^pi^qi ^Pj^qj> 

p,q = 1, m. Similarly we define n x n matrix op- 
erators ro,Ti,r2 and 1^ by replacing the indices i,j 
in Eo,Ei,E2 and E3 with k, I respectively, and setting 
p, q = 1, n. It is straightforward to derive that C(p a p) 
can be expressed as the mean values of the above local 
observables, 



C 2 ( Pa p) = 



i + ^((E3®r 3 ) 2 -(S3®r ) 2 



(S C 
(So- 



r 2 ) 2 



(S ( 

(s 3 



Ti> 2 - 

>r 2 ) 2 ) 



Tr) 2 
(16) 



In summary we have presented an experimental way 
to measure the entanglement of formation for arbitrary 
dimensional pure states, by measuring some local quan- 
tum mechanical observables. We reduced the difficult 
problem to find the concurrence of "two-qubit" states 
for which many results have been already derived. Re- 
cently high dimensional bipartite systems like in NMR 
and nitrogen-vacancy defect center have been success- 
fully used for quantum computation and simulation ex- 
periments Our results present a plausible way to 
measure the entanglement of formation in these systems 
and to investigate the roles played by the entanglement 
of formation in these quantum information processing. 

So far experimental measurement on entanglement of 
formation and concurrence concerns only pure states. For 
mixed states, less is known except for experimental deter- 
mination of separability, both sufficiently and necessary, 
for two-qubit and qubit-qutrit systems [2(| • Gener- 
ally (j4]) has only analytical results for some special states 
[2l| and analytical lower bounds |23| which are not exper- 
imentally measurable. Recently in [23[ we have presented 
a measurable lower bound of entanglement of formation. 
The formula (HJ for pure state may also help to study 
measurable lower bounds of entanglement of formation 
for mixed states. 
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